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ADDITIVE  AND  MULTIPLICATIVE  DISTORTION 
SIMULATOR  FOR  NON-REAL  TIME  TESTING 
OF  HF  MODEMS 

1 .0  Introduction 

This  report  describes  an  hf  channel  simulator  which  includes  the  propagation 
path,  the  frequency  response  of  the  rf  equipment,  inband  interference  from  other 
users  of  a  common  channel,  and  the  operational  conditions.  The  operational 
conditions  are  related  primarily  to  doppler  problems  incurred  with  high  velocity 
aircraft. 

This  is  an  off-line,  non-real  time  test  capability  using  sampled  data 
input  from  a  modulator  and  providing  sampled  data  output  to  a  demodulator. 

Test  conditions  are  established  in  a  table  of  variables,  permitting  a  flexible 
means  of  selecting  combinations  of  conditions.  It  is  the  purpose  of  this  report 
to  describe  the  capabilities  of  this  simulation  program. 

2.0  Simulation  Program 
2.1  Types  of  Stress 

A  FORTRAN  program  (Appendix  A)  has  been  developed  which  permits  the 
application  of  a  set  of  selectable  stress  conditions  to  a  digitized  modem  line 
signal.  The  program  is  independent  of  the  modem  signal  design,  depending  only 
on  the  sampling  rate  used  to  digitize  the  input  signal. 

The  types  of  stress  for  which  test  conditions  may  be  specified 

include: 

.  Gaussian  noise  level 

.  frequency  translation  error 

.  multipath  delay  with  fixed  doppler  offset 

/ 

.  varying  doppler  due  to  operational  scenario 
.  multi  tone  interference 

Manuscript  submitted  April  9,  1981. 
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.  random  noise  pulse  interference 
.  equipment  filter  characteristics 
.  AGC  characteristics. 

The  program  uses  independent  subroutines  for  introducing  each  stress 
condition.  The  parameters  for  describing  the  test  conditions  are  contained  in 
a  two-dimensional  (10  x  10)  floating  point  array  which  is  accessed  by  the  sub¬ 
routines.  Table  1  is  a  list  of  all  the  variables  in  that  array  with  a  description 
of  their  use.  Appendix  B  is  a  FORTRAN  program  used  to  generate  the  table,  and 
Appendix  C  is  a  program  to  modify  an  existing  table. 

Figure  1  is  a  block  diagram  showing  the  order  in  which  the  stress 
functions  may  be  performed.  Any  of  the  functions  may  be  readily  bypassed.  Noise 
and  interference  are  added  after  introducing  the  multipath  and  doppler  distortions. 
Filtering  is  divided  into  a  transmit  and  a  receive  filter,  making  it  possible  to 
filter  the  desired  signal  twice,  but  the  interference  only  once.  Table  2  is  a 
list  of  the  subroutines  and  their  use. 

2.2  Input/Output 

The  test  program  operates  from  a  sampled  digitized  input  of  the  modulator 
output  signal.  The  input  may  be  from  either  magnetic  tape  or  a  disk  file.  In 
either  case,  each  sample  is  quantized  to  12  bits  and  left  justified  in  a  16  bit 
word.  The  data  are  stored  in  blocks  of  1000  samples.  The  magnetic  tape  is  a 
standard  nine-track,  800  byte  per  inch  recording  system. 

Each  integer  input  sample  is  converted  to  floating  point  and  normalized 
to  a  value  corresponding  to  zero  dBm  across  600  ohms.  Normalization  of  the  input 
data  permits  all  interference  levels  to  be  specified  relative  to  zero  dBm. 

The  sampled  data  is  sequenced  through  the  various  subroutines  in  blocks 
of  equal  size,  whose  length  is  specified  in  the  table  of  variables.  It  is 
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Table  1  —  Identity  of  variables  in  control  table 
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Fig.  1  —  Non-real-time  stress  program 
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Table  2  —  Identity  of  subroutines 


1. 


2. 


3. 


4. 


5. 


6. 


7. 


8. 


9. 


10. 

11. 

12. 

13. 

14. 

15. 


INPUTX 

Read  in  NBD  sample  of  data  from  either  maonetic  tape  or  a  disk  file 
representing  modulator  output  signal.  ' 

Tx  FILT 

Transmitter  filter;  selection  of  one  of  five  algorithms  to  filter 
input  data. 

ALPASX 

All -pass  filter  NBD  samples  of  data  to  generate  inphase  and  quadrature 
components. 

DELAYX 

Combine  output  of  N  propagation  paths  for  different  delays. 

FREQX 

Frequency  translate  composite  signal  by  a  fixed  or  variable  amount. 
OOPVAR 

Update  doppler  frequency  in  accordance  with  one  of  three  operational 
scenarios. 

INTERF 

Add  non-fading  tonal  interference  to  fading  signal.  Interference  may 
be  CW,  FSK,  DPSK,  and/or  Swept  Tone. 

PULSE 

Add  wideband  Gaussian  noise  pulse  to  received  signal. 

NOISEX 

Add  Gaussian  noise  to  received  signal. 

RxFILT 

Receiver  filter,  selection  of  one  of  five  algorithms  to  filter  output 
signal  plus  noise. 

AGC 

Receiver  automatic  gain  control. 

STOREX 

Store  NBD  samples  of  data  in  integer  record  for  transfer  to  either 
magnetic  tape  or  disk  file. 

MTIN 

Read  block  of  1000  samples  from  magnetic  tape. 

MTSET 

Set  MTOUT  for  D/A  operation 
MTOUT 

Write  block  of  1000  samples  on  magnetic  tape. 
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convenient  to  make  this  block  size  equal  to  the  number  of  samples  in  a  baud  of 
the  modulated  tonal  interference.  In  the  absence  of  such  an  interference,  the 
block  size  may  be  any  value  up  to  100  samples. 

After  all  signal  processing,  the  normal ization  is  removed  and  the 
distorted  signal  is  converted  to  integer  and  stored  on  the  output  medium  in  the 
same  format  as  the  input  signal.  The  output  medium  may  be  either  magnetic  tape 
or  a  disk  file. 


2.3  Transmission  Equipment  Filters 


Provisions  are  made  to  select  one  of  five  forms  of  a  simple  non¬ 
recursive  bandpass  filter.  Characteristics  of  the  five  filters  are  completely 
specified  within  the  filter  subroutines.  All  five  filters  are  symmetrical  about 


the  mid-band  frequency,  which  is  one- fourth  of  the  sampling  rate.  The  filters 
are  defined  by  the  following  equations: 


1 . 

y0 

-  0.5  X0  -  0.5  X_2 

2. 

y0 

=  0.25  XQ  -  0.5  x_2 

+  0.25  X  , 

3. 

^0 

=  0.6  X0  -  0.5  X.2  - 

■  0.1  X_4 

4. 

y0 

=  0.36  X0  -  0.6  X_2 

+  0.13  X__4  +  0.1  X 

5. 

y0 

=  0.1296  X0  -  0.432 

X_2  +  0.453  X.4  - 

-0.0959  X_8  +  0.014  X.io  +  0.0127  X_i2 
+0.002  X_14  +  0.0001  X_16 

Table  3  gives  the  amplitude  and  phase  response  of  the  five  filters. 


Filter  one  is  the  simplest  form  of  a  non-recursive  bandpass  filter.  It  has  a 
3  dB  bandwidth  of  0.5  W,  where  W  is  the  Nyquist  frequency.  The  phase  response 
is  linear.  Filter  two  is  a  cascade  of  two  of  the  number  one  filters.  The  result 
is  a  6  dB  bandwidth  of  0.5  W  and  a  3  dB  bandwidth  of  0.36  W.  The  phase  response 
is  still  linear. 


6 


Table  3  —  Amplitude  and  phase  response  of  bandpass  filters 


FMter  =1 

Filter  *2 

Filter  43 

Filter 

:iA 

Filter  «5 

AMP 

PHASE 

AMP 

PHASE 

AMP 

PHASE 

AMP 

PHASE 

AMP 

PHASE 

Kl 

.0000 

90.00 

.0000 

90.00 

.0000 

.0005 

SH 

.1564 

81.00 

.0245 

162.00 

.2177 

78.45 

.0474 

156.91 

.0017 

-52.69 

.3090 

72.00 

.0955 

144.00 

.4224 

.1784 

134.13 

.0312 

-92.26 

.15 

.4540 

63.00 

.2061 

126.00 

.6027 

KtKSiS 

.3632 

.1312 

-136.12 

.20 

.5878 

54.00 

.3455 

108.00 

.7500 

45.43 

.5626 

90.85 

.3160 

478.33 

.25 

.7071 

45.00 

.5000 

90.00 

.3602 

35.54 

.7400 

.5472 

142.12 

.30 

.3090 

36.00 

.6545 

72.00 

.9336 

26.51 

.8716 

53.03 

.7593 

106.02 

.35 

.3910 

27.00 

.7939 

54.30 

.9752 

13.50 

.9510 

37.00 

.  9033 

73.97 

.40 

.9511 

18.00 

.9045 

36.00 

.9937 

11.54 

.9874 

WmiI 

.9744 

46.14 

.45 

.9377 

9.00 

.9755 

18.00 

.9992 

5.50 

VH6t1 

.9963 

21.99 

.50 

1.0000 

0.00 

1.0000 

0.00 

1.0000 

0.00 

1 . 0000 

0.00 

.3995 

0.00 

.55 

.9377 

-9.00 

.9755 

-18.00 

.9992 

-5.50 

.9984 

.9963 

-21.99 

.60 

,9511 

-18.00 

.9045 

-36.00 

.9937 

-11.54 

.9874 

-23.08 

.9744 

-46.14 

.65 

.3910 

-27.00 

.7939 

-54.00 

.9752 

-18.50 

.9510 

-37.00 

.  3033 

-73.97 

.70 

.3090 

-36.00 

.6545 

-72.00 

.9336 

-26.51 

.8716 

-53.03 

.7593 

-106.02 

.75 

.7071 

-45.00 

.5000 

-90.00 

.3602 

-35.54 

-71.08 

.5472 

-142.12 

.30 

.5373 

-54.00 

. 3455  • 

■108.00 

.7500 

-45.43 

.5626 

-90.35 

.2160 

178.33 

.35 

.4540 

-53.00 

.2061 

■126.00 

.6027 

SriTiiH 

.3632 

-112.00 

.1312 

136.12 

.90 

.3090 

-72.00 

.0955  ■ 

•144.00 

.4227 

RtXrfl 

.1784 

-134.13 

.0312 

92.26 

.95 

.1564 

-81.00 

.0245  ■ 

•162.00 

.2177 

-78.45 

.0474 

-156.9! 

.0017 

52.69 

1.00 

.0000 

-90.00 

.0000 

-90.00 

.0000 

-79.65 

-180.00 

.0005 

-180.00 

Filter  three  has  a  flat  top  amplitude  response  with  a  3  dB  bandwidth 
of  0.63  W.  The  phase  response  is  still  quite  linear.  This  filter  is  equivalent 
to  a  parallel  implementation  of  two  filters,  yQ  =  0.5  x0  -  0.5  X  2  and  yQ  = 
o.i  XQ  -  o.l  X_4 

Filter  four  is  a  cascade  of  two  of  the  number  three  filters.  The 
result  is  a  6  dB  bandwidth  of  0.63  W  and  a  3  dB  bandwidth  of  0.51  W.  Filter 
five  is  a  cascade  of  four  of  the  number  three  filters.  The  result  is  a  12  dB 
bandwidth  of  0.63  W,  a  6  dB  bandwidth  of  0.51  W,  and  a  3  dB  bandwidth  of  0.43  W. 

For  convenience,  a  separate  filter  subroutine  is  designated  for  the 
transmitter  and  the  receiver  filter  functions.  Two  values  in  the  table  of 
variables  control  which  filter  algorithms  are  selected.  These  values  may  be 
zero  through  five,  with  zero  representing  no  filtering.  It  is  not  necessary 
that  the  same  filter  characteristics  be  used  for  both  the  transmit  and  receive 
fil ters . 

2.4  Analytic  Signal 

The  envelope  of  the  real  signal  is  the  absolute  value  of  the  signal 
preenvelope  (analytic  signal), 
fa(t)  =  f(t)  +  jf(t) 

where 

fa(t)  is  the  analytic  signal  as  a  function  of  time 

f(t)  is  the  real  signal 

f(t)  is  the  Hilbert  transform  of  f(t) 

The  benefit  of  the  analytic  signal  is  that  it  provides  a  means  to  frequency 
shift  the  transmitted  signal  in  the  time 'domain.  This  may  be  performed  by 
multiplying  the  analytic  signal  by  a  complex  exponential. 
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A  digital  all-pass  filter  is  used  to  generate  an  approximation  to  the 
ideal  Hilbert  transform  pair.  The  input  to  the  all-pass  filter  is  the  real  signal. 
The  two  outputs  are  the  in-phase  and  quadrature  components  of  the  analytic  signal, 
that  have  nearly  equal  magnitude  and  are  approximately  90  degrees  apart.  An 
alternate  method  of  generating  an  approximation  to  the  ideal  Hilbert  transform 
is  by  means  of  a  non-recursive  discrete  Hilbert  transform.  Such  designs  can 
yield  an  ideal  phase  response,  but  there  are  fluctuations  in  the  magnitude 
response,  depending  on  the  number  of  terms  used. 

2.5  Multipath 

Program  inputs  from  the  table  of  variables  permit  the  user  to  specify 


the  number  of  propagation  paths,  their  relative  delays,  and  their  relative 
amplitudes.  The  mean  frequency  offset  and  the  doppler  spread  of  each  path  are 
determined  directly  from  an  input  data  file.  One  means  of  generating  that  data 
is  by  the  FORTRAN  program  in  Appendix  D.  That  program  generates  the  complex 
function  representing  the  fading  on  each  path  at  a  sampling  rate  equal  to  1/100 
of  the  modem  sampling  rate.  Each  value  is  equal  to  the  i nstantaneous  amplitude  of  the 
sum  of  several  sinusoids  with  small  differences  in  frequency. 


E(t)  =  ^  cos(w  +  wn)  t  +  j  sin  (w  +  wn; 
n=l 


where 


w  is  the  mean  frequency  offset  of  the  path;wn  is  one  of  N  frequencies 
selected  to  be  approximately  (but  not  exactly)  uniformly  distributed 
over  the  fading  bandwidth 

The  values  for  each  path  are  interleaved  and  stored  in  blocks  of  1024  samples. 
These  values  are  generated  at  1/100  of  the  modem  sampling  rate  in  order  to 
lessen  the  storage  requirements.  This  channel  sampling  rate  is  well  above  the 
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minimum  rate  requirements  for  the  fading  channel.  The  simulation  program  inter¬ 
polates  between  input  values  to  obtain  values  at  the  modem  sampling  rate.  Linear 
interpolation  is  used. 

The  program  in  Appendix  D  for  generating  sampled  values  of  the  fading 
on  the  propagation  paths  uses  N  equal  amplitude  complex  sinusoids,  where  N  is 
a  variable.  The  frequency  of  each  sinusoid  is  chosen  randomly  with  the 
restriction  that  there  is  one  assignment  in  each  of  N  equal  width  increments 
over  the  total  range  specified  for  the  frequency  spread  due  to  the  propagation 
path.  The  inputs  to  that  program  are  the  number  of  tones  per  paths,  the 
number  of  paths,  the  frequency  spread  for  each  path,  and  the  mean  frequency 
offset  for  each  path.  The  relative  amplitude  of  the  different  paths  is  a 
variable  in  the  main  simulation  program. 

The  complex  gain  values  are  generated  by  combining  N  sinusoids,  each 
of  which  have  a  rms  amplitude  of  1.0.  Thus,  the  rm$  amplitude  of  the  composite 
signal  is  equal  to  fTT.  To  keep  the  rms  amplitude  of  the  composite  signal  inde¬ 
pendent  of  the  number  of  sinusoids  used  to  generate  it,  the  complex  components 
are  divided  by  i^N.  Thus,  the  normalized  rms  amplitude  of  the  composite  signal 
is  always  1.0. 

Figure  2  is  a  plot  of  the  cumulative  distribution  of  the  instantaneous 
amplitude  of  the  fading  signal  for  a  one-path  channel  model  with  a  spread  of 
1.0  Hz.  These  data  were  generated  by  combining  20  sinusoids.  The  frequencies 
of  the  20  sinusoids  are  listed  in  Table  4.  The  modem  sampling  rate  was  specified 
as  7200  Hz,  thus  the  sample  rate  for  the  path  data  was  72  Hz.  Figure  3  is  a 
plot  of  the  median  fade  duration  vs.  doppler  spread.  A  fade 

was  defined  as  the  period  of  time  for  which  the  instantaneous  amplitude  was 
below  a  threshold.  More  exactly,  the  beginning  of  a  fade  was  defined  as  when 
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Table  4  —  Example  of  random  tone  selection  for  spread  of  1  Hz 


NO. 

FREQ.  (Hz) 

NO. 

FREQ.  (Hz) 

1 

'0. 951 

11 

0.079 

'  2 

-a  878 

12 

0.169 

3 

-0.790 

13 

0.277 

4 

-a  638 

14 

a  381 

5 

-a  594 

15 

a  420 

6 

-0.404 

16 

0  599 

7 

-a  353 

17 

0.633 

8 

-0. 261 

18 

0.763 

9 

-0.129 

19 

0  802 

10 

-0.048 

20 

0.994 

12 


MEDIAN  FADE  DURATION  (SEC) 


the  level  of  the  signal  decreased  to  -6  dB  relative  to  the  mean,  and  the  end 
of  a  fade  was  defined  as  when  the  level  of  the  signal  increased  to  -5  dB 
relative  to  the  mean.  Fade  rate  was  defined  as  the  average  number  of  times 
the  signal  level  left  a  fade  condition  in  one  minute.  For  the  channel  model 
described  above,  the  average  fade  rate  was  45  fades/minute  measured  over  a  10 
minute  period.  Figure  4  is  a  plot  of  fade  rate  versus  spread  for  a  one-path 
model  combining  20  sinusoids  at  a  modem  sampling  rate  of  7200  Hz. 

The  simulation  program  provides  for  specifying  up  to  ten  paths.  The 
path  delays  are  specified  in  seconds,  but  they  are  interpreted  in  increments  of 
the  sampling  period.  The  maximum  amount  of  delay  cannot  exceed  the  period  of 
the  block  of  data  processed  by  each  pass  through  the  subroutines.  Thus,  if  data 
were  processed  in  blocks  of  100  samples  at  7200  samples/second,  the  maximum 
delay  would  be  0.01389  seconds. 

The  table  of  variables  contains  the  relative  rms  amplitude  of  each 
propagation  path.  The  program  calculates  an  amplitude  scale  factor  for  each 
propagation  path  which  results  in  the  received  composite  signal  having  the  same 
rms  level  as  the  transmitted  signal.  The  scale  factor  for  each  path  is  computed 
as 


where 

Ei  is  the  rms  signal  level  for  signal  i 
2.6  Popp! er 

In  addition  to  the  fixed  doppler  shift  for  each  propagation  path,  the 


program  provides  for  a  variable  doppler  representing  the  composite  of  the  frequency 


translation  errors  and  doppler  shift  due  to  the  relative  velocity  between  the 
transmitter  and  receiver  platforms.  This  composite  frequency  error  may  be 
either  fixed  or  time  varying.  Provisions  are  made  for  the  selection  of  one  of 
three  types  of  varying  doppler.  They  are: 

A  step  change  in  doppler  of  a  specified  amount  to  occur  at 
a  specified  time  after  the  start  of  the  program.  Time  is 
computed  in  increments  of  the  sampling  period.  This  vjye 
of  test  represents  the  condition  where  there  is  an  extremely 
rapid  change  in  the  relative  velocity  between  the  transmitter 
and  receiver  after  a  communication  link  has  been  established. 
It  also  may  represent  a  step  change  in  the  frequency  controls 
of  the  radio  equipment.  It  provides  a  means  to  measure  the 
ability  of  the  modem  receiver  to  accommodate  a  step  change 
in  frequency. 

A  ramp  change  in  doppler  at  a  specified  rate  with  the 
direction  of  the  ramp  reversing  when  high  and  low  limits 
are  reached.  Each  incremental  change  in  doppler  is  equal 
to  the  rate  of  change  divided  by  the  sampling  rate.  The 
linearly  changing  doppler  with  time  provides  a  means  to  test 
the  doppler  tracking  capability  of  the  modem  receiver.  Also, 
by  changing  the  high  and  the  low  limits,  it  provides  a  means 
to  test  the  effects  of  variations  in  the  amplitude  and  phase 
delay  due  to  the  response  of  the  receiver  filter  at  the  band 
edges . 

A  high  velocity  fly-by  scenario  as  diagramed  in  Figure  5. 

For  this  test  the  specifications  include  the  velocity  of 


0.5 


1.0 


2.0 


4.0 


SPREAD  (Hz) 

Fade  rate  vs.  frequency  spread,  one  path,  software  channel  and  CSP-30 
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1 


Fig.  5  —  Fly-by  scenario 


the  aircraft,  the  radio  frequency,  the  minimum  distance 


reached  between  transmitter  and  receiver  at  time  tQ,  the 

time  (ti)  before  t0  at  which  to  start  computing  the  doppler 
change,  and  that  time  relative  to  the  start  of  the  program 
(e.g.,  three  seconds  after  the  start  of  the  program, 
initialize  the  fly-by  test.  At  that  time,  the  aircraft 
will  be  ten  seconds  from  tQ.  At  t0  the  distance  between  the 
transmitter  and  the  receiver  will  be  two  miles.  The  aircraft 
velocity  is  1000  miles  per  hour  and  the  radio  frequency  is 
10  MHz). 

The  general  equation  is: 

Doppler  =  cos  (6) 

=  i£  y 

C  (x2  +  y2)* 

=  V2F 

C  (xz  +  (V(tx  -t))2)*s 

With  the  aircraft  velocity  expressed  in  miles/hour,  the  radio  frequency 
in  MHz,  the  distance  in  miles,  and  time  in  seconds,  the  doppler  may  be 
calculated  as  follows: 

Doppler  =  ( — L.)2  _ F  _ 

3600  0.18b  ^x2+(v(t  -t))2^ 

3600 

where : 


Doppler  is  in  Hertz 
V  is  velocity  in  mph 
F  is  frequency  in  MHz 

t^  is  time  away  from  tQ  at  start  of  test,  in  seconds 
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t  is  elapsed  time  in  test,  in  seconds 

X  is  distance  between  transmitter  and  receiver  at  time  tQ,  in  miles. 

Table  5  is  a  tabulation  of  the  resulting  doppler  at  increments  of  10 
seconds  for  the  fly-by  conditions  given  in  the  previous  example.  For  this 
scenario,  the  maximum  stress  on  the  modem  doppler  tracking  capability  is  at  time 
t0.  Figure  6  shows  the  change  in  doppler  at  time  tQ  versus  aircraft  velocity  and 
distance  between  transmitter  and  receiver  for  a  propagation  frequency  of  10  MHz. 

2.7  Tonal  Interference 

Interference  corresponding  to  other  users  of  the  radio  channel  is 
provided  for  by  specifying  tonal  interference  to  be  either  CW,  binary  FSK,  four- 
phase  DPSK,  or  a  saw-tooth  swept  audio  tone.  The  number  of  tones,  their  frequency 
assignments,  modulation  baud  periods,  and  rms  power  level  are  specified.  If  the 
tones  are  modulated,  the  information  is  random  data,  but  the  baud  period  must 
correspond  to  the  same  period  as  the  block  of  data  being  processed. 

The  interference  may  be  any  combination  of  the  four  types.  If  two 
different  types  of  modulated  interference  are  specified,  the  conditions  are 
limited  to  that  in  which  they  are  both  simultaneously  keyed  at  an  identical 
rate. 

The  interference  is  added  to  the  signal  after  summation  of  the  signals 
from  each  propagation  path,  and  after  frequency  shifting  all  components.  Thus, 
the  interference  is  neither  distorted  by  multipath  nor  is  it  frequency  shifted. 

This  represents  the  condition  where  there  is  always  a  stationary  path  between  the 
interfering  transmitter  and  modem  receiver.  For  the  case  of  a  high  velocity 
fly-by  scenario  with  tonal  interference,  it  represents  transmission  of  the  modem 
signal  from  the  aircraft  to  the  ground  station,  because  receiving  by  the  aircraft 
would  require  a  doppler  shift  of  the  interference.  Likewise,  if  there  is  a 
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Table  5  —  Doppler  Incurred  in  fly-by  example,  Vel  =  1000  mph, 
RF  =  10  MHz,  minimum  distance  =  2  miles 


TIME 

BEFORE 

,To  . 

(sec) 

DOPPLER 

(Hz) 

DOPPLER 

CHANGE 

(Hz/sec) 

10 

7.36 

9 

7.25 

-0.11 

8 

7.11 

-0.14 

7 

6.92 

-0.19 

6 

6.65 

-0.27 

5 

6.27 

-0.38 

4 

5.71 

-0.56 

3 

4,88 

-0.83 

2 

3.67 

-1.21 

1 

2.01 

-1.66 

0 

0.00 

-2.01 

0  500  1000  1500  2000  2500  3000 

VELOCITY  (mph) 

Fig.  6  —  Doppler  change  during  last  second  during  fly-by  scenario,  as  a  function  of  velocity 
and  distance,  for  a  transmission  frequency  of  10  MHz 
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fixed  frequency  translation  error  plus  tonal  interference,  the  translation  error 
must  be  due  to  the  signal  transmitter.  If  it  is  desired  to  represent  a  condition 
where  the  interfering  transmitter  has  a  frequency  translation  error,  this  may  be 
done  by  making  a  corresponding  change  in  the  frequency  assignments  for  the 
i nterference . 

The  level  of  the  tonal  interference  is  specified  in  dBm  across  600  ohms. 
The  program  calculates  the  root  mean  square  (RMS)  voltage  of  the  composite  signal 
as  follows: 

as,  V(dBm)  =  20  log  (V(RMS )/ yjoTT) 

then,  V(RMS)  =  .JO-  10  v(dBm)/20 

This  is  normalized  to  the  RMS  voltage  for  zero  dBm  so, 

V1  (RMS)  =  V(RMS)  /yjOT~ 

=  10V(dBm)/20 

The  peak  amplitude  of  each  interference  tone  is  related  to  the  RMS  amplitude  of  the 
composite  signal  as  follows: 

E(peak)  =  V'  (RMS)/^72 

=  V2/N  10  V(dBtT!)/20 

The  four  types  of  tonal  interference  may  occur  in  any  combination  as 
specified  by  an  entry  in  the  table  of  variables.  The  program  converts  that  value 
from  a  floating  point  number  to  an  integer  with  values  of  zero  through  15.  Zero 
indicates  no  interference  is  to  be  generated.  The  four  least  significant  bits 
of  that  integer  control  the  types  of  interference  to  be  generated  as  follows: 

Bit  No.  Set  to  1  Interference  Generated 

0  CW 

1  FSK 

2 


3 
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DPSK 

Swept  Tone 


Table  6  lists  the  control  word  for  each  of  the  15  possible  combinations  of  tonal 
interference . 

2.8  Pulse  Interference 

Provisions  are  made  in  the  simulation  program  for  a  subroutine  for  the 
addition  of  non-Gaussian  noise.  At  the  present  time  an  acceptable  routine  has 
not  been  developed.  A  dummy  subroutine  (Subroutine  Pulse)  is  included  with 
provisions  for  characterizing  the  noise  pulses  in  terms  of  their  RMS  power 
level,  duration,  and  the  probability  of  occurence. 

2.9  Gaussian  Noise 

The  density  of  the  continuous  additive  Gaussian  noise  is  specified  in 
dBm/Hz.  The  program  converts  that  value  to  a  noise  scale  factor  which  is  used  to 
multiply  each  Gaussian  noise  sample.  Each  noise  sample  is  generated  from  a  pair 
of  random  numbers,  with  uniform  distribution  between  zero  and  one.  The  equation  is 

X  =  ANZ  (-2.0  1  n(U ) )  cos  (2ttV) 

where : 

X  =  a  sample  of  Gaussian  noise 
ANZ  =  noise  scale  factor 

U  and  V  =  two  uniformly  distributed  random  numbers  (zero  to  one). 

With  a  noise  scale  factor  of  1.0,  the  algorithm  generates  Gaussian 
noise  with  a  long-term  RMS  value  of  1.0  (zero  mean  and  variance  =1).  The 
scale  factor  required  to  generate  the  specified  noise  density  level  is  calculated 
as  follows: 

P(dBm/Hz )  =  20  log  (ANZ)  -  10  1 og (SAMP/2 ) 
where:  P(dBm/Hz)  is  the  specified  noise  density 

ANZ  is  the  noise  scale  factor 

SAMP  is  the  modem  sampling  rate 

( P(dBm/Hz )  +  10  log  (SAMP/2  )  )/20 

s0’  ANZ  =  10 
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Table  6  —  Control  words  for  selection  of  tonal  interference 


CONTROL 

WORD 

(OCTAL) 

TYPES  OF  INTERFERENCE 

00 

None 

01 

CW 

02 

FSK 

03 

CW,  FSK 

04 

DPSK 

05 

CW,  DPSK 

06 

FSK,  DPSK 

07 

CW,  FSK,  DPSK 

10 

SWEPT  TONE 

11 

CW,  SWEPT  TONE 

12 

FSK,  SWEPT  TONE 

13 

CW,  FSK,  SWEPT  TONE 

14 

DPSK,  SWEPT  TONE 

15 

CW,  DPSK,  SWEPT  TONE 

16 

FSK,  DPSK,  SWEPT  TONE 

17 

CW,  FSK,  DPSK,  SWEPT  TONE 

2.10  AGC 

The  radio  receiver  automatic  gain  control  characteristics  may  have  a 
significant  effect  on  the  performance  of  some  modems.  Provisions  are  made  in 
the  simulation  program  for  an  AGC  function,  but  at  the  present  time  this  is  just 
a  dummy  subroutine.  The  table  of  variables  include  the  following  AGC  variables: 
Var  (2,2)  -  current  AGC  gain  value 
Var  (5,2)  -  constant  controlling  AGC  attack  time 
Var  (6,2)  -  constant  controlling  AGC  decay  time 
Var  (7,2)  -  constant  representing  maximum  AGC  amplification 

2.11  Input  Variables  Modified  by  Program 

Each  test  is  controlled  by  the  variables  read  from  a  disk  file  into 
the  two  dimensional  array  VAR  (10,  10).  These  values  are  all  entered  as  floating 
point  numbers.  The  following  is  a  list  of  those  variables  that  are  converted  to 
integer  by  the  program,  but  otherwise  remain  unchanged  during  the  entire  test: 

1 .  VAR  (2,  1 )  =  N8D 

2.  VAR  (3,  1)  =  NPATH 

3.  VAR  (4,  1)  =  NFADE 

4.  VAR  (5,  1)  =  KDOPP 

5.  VAR  (6,  1)  =  KINTF 

6.  VAR  (7,  1)  =  KPULSE 

7.  VAR  (8,  1)  =  KTXFIL 

8.  VAR  (9,  1)  =  KRXFIL 

9.  VAR  (10,  1)  =  KAGC 

10.  VAR  (1,5)=  NCW 

11 .  VAR  (1,6)=  NFSK 

12.  VAR  (1,7)=  NPSK 

13.  VAR  (3,  2)  =  NTAPE 

14.  VAR  (4,  2)  =  MTAPE 
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The  interference  levels  of  the  tonal  interference  are  specified  in 
dBm.  The  program  converts  these  to  the  peak  amplitude  of  each  tone.  The 
variables  involved  are: 

1 .  VAR  (4,  5)  =  AMPCW 

2.  VAR  (5,  6)  =  AMPFSK 

3.  VAR  (4,  7)  =  AMPPSK 

4.  VAR  (4,  8)  =  AMPSAW 

The  interference  levels  for  the  Gaussian  noise  and  the  noise  pulse  are 
specified  in  dBm/HZ.  The  program  converts  these  values  to  scale  factors.  The 
variables  involved  are: 

1 .  VAR  (1,2)=  GAUSS 

2.  VAR  (2,  10)  =  AMPPUL 

The  path  delays  are  specified  in  seconds.  The  program  multiplies 
these  by  the  sampling  rate  to  get  the  delay  in  terms  of  the  number  of  sample 
periods.  The  variables  involved  are  in  VAR(I,  3),  where  I  refers  to  the  path 
number. 

The  propagation  paths  are  defined  by  their  relative  RMS  amplitudes. 

The  program  converts  these  to  scale  factors  dependent  on  the  number  of  paths. 

The  variables  affected  are  in  VAR  (I,  4),  where  I  refers  to  the  path  number. 

3.0  Conclusions 

A  FORTRAN  program  has  been  developed  to  provide  a  simple,  standardized 
means  of  applying  uniform  tests  to  various  hf  modem  signals.  It  has  been  applied 
to  testing  of  the  ANDVT  hf  modem  signal  design,  the  results  of  which  are  published 
in  reference  [16].  The  program  provides  a  beneficial  supplement  to  real-time 
testing  and  to  Monte-Carlo  type  non-real -time  ti.-i.-s. 
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APPENDIX  A 
PROGRAM  STRESS. FTN 
TO  ADD  STRESS  TO  MODEM  SIGNAL 
IN  FORM  OF  GAUSSIAN  NOISE,  DOPPLER,  FILTERING 
NARROWBAND  INTERFERENCE,  PULSE,  MULTIPATH 

INPUT  IS  DIGITIZED  MODULATOR  OUTPUT  ON  MT:  OR  DISK 
FILE  TAN. DAT 

OUTPUT  IS  DIGITIZED  INPUT  TO  DEMODULATOR 
ON  MT:  OR  DISK  FILE  STRESS . DAT 

REQUIRED  DATA  FILES  FROM  DISK  ARE: 

STRVAR.DAT  TABLE  OF  VARIABLES 
PATH . DAT  RAYLIEGH  PATH  DATA 

COMMON /BLKl/VARC 10, 10 ) , A( 1 00 ) , Q( 1 00 ) 

COMMON /BLK2/TXFIL (200) 

COMMON/BLK3/DELAY(2, 200) 

C0MM0N/BLK4/PHASE (40 ) , TAU ( 4  )  , FQSAW 
COMMON/BLK5/PATH(1024) ,WT, TIMEX 
C0MM0N/BLK6/I 1 , 12 
C0MM0N/BLK7/FLAG, KFLAG 
C0MM0N/BLK8/IA( 1 000 ) , ISTAT , NDX , BLOCK 
C0MM0N/BLK9/SAVIN ( 3,2), SAVOUT (3,2) 

COMMON/BLK10/RXFIL(200) 

C0MM0N/BLK1 1/AGCFIL(2, 100) 

C0MM0N/BLK12/IB( 1000) , JSTAT , NDZ , KCODE , RMS , PE AK , TOTAL 
C0MM0N/BLK1 3/DELTA, ROLD(10 ) ,RNEW( 10) ,GOLD( 10) ,GNEW( 10) ,NDY,KSTAT 
COMPLEX  PATH 

EQUIVALENCE  (VAR(1 , 1 ) , SAMP ) , ( VAR (2 , 1 0 ) , AMPPUL ) , ( VAR ( 1 ,2), GAUSS) 
EQUIVALENCE  (VAR (4 , 5 ) , AMPCW)  , ( VAR ( 5 , 6 ) , AMPFSK )  , ( VAR ( 4 , 7 ) , AMPPSK ) 
EQUIVALENCE  (V AR ( 1 ,8),F1SAW) , (VAR ( 4 , 8 ) , AMPSAW) 

WRITE (5 , 1 ) 

1  FORMAT( IX, 'ENTER  NO.  OF  FRAMES',/) 

READ(5 , 2  )FRAMES 

2  FORMAT (F8 . 0 ) 

WRITE (5,3) 

3  FORMATOX, 'ENTER  NO.  OF  BLOCKS  TO  SKIP,  STARTING  POINT',/) 

READ(5 , 4 )NSK 
READ(5 , 4  )NDX 

4  FORMAT (15) 

NSK=NSK+1 

C  READ  TABLE  OF  VARIABLES  FROM  DISK 

CALL  ASSIGN(4, 'STRVAR.DAT') 

READ(4,END  =  3000, ERR =3000  )VAR 
CALL  CLOSE  (4) 

C 

C  SET  INITIAL  CONDITIONS 

NBD=VAR (2 ,  1  ) 

NPATH=VAR (3  ,  1  ) 

NFADE  =VAR ( 4 ,  1  ) 


A-l 


o  o 


5 


c 


7 

C 

C 


C 


9 


1 1 

C 

c 

15 


KDOPP=VAR (5  ,  1  ) 

KINTF  =VAR (6 , 1  ) 

KPULSE=VAR (7  ,  1  ) 

KTXF IL=VAR (8 ,  1  ) 

KRXFIL=VAR(9, 1  ) 

KAGC  =VAR ( 1 0 , 1  ) 

NCW=VAR ( 1 ,5) 

NFSK=VAR ( 1 ,6) 

N  PSK=V  AR ( 1 ,7) 

NTAPE  =VAR ( 3 , 2 ) 

MTAPE  =VAR (4 , 2 ) 

DO  5  1=1,4 
TA U (I ) =0  . 

KFLAG=0 
RMS=0 . 

PEAK=0. 

T0TAL=0 . 

FLAG=0 . 

11=0 
12  =  0 

FRAME=0 . 

ND2  =0 

FQSAW=F 1SAW 
WT=0  . 

KC0DE=3 
TIMEX  =0 . 

RANDOM  PHASES  FOR  DPSK  TONES 

DO  7  1=1 , 40 
X=RAN(I1 , 12) 

PHASE(I)  =  2.0*(X-0.5)«3.  14159 

ASSIGN  AND  SET  INPUT  DEVICE 
NTAPE=4  IS  NO  INPUT 

IF(NTAPE.EQ.4)G0  TO  15 
IF(NTAPE.EQ.3)G0  TO  10 

INPUT  FROM  MT: 

CALL  ASSIGN ( 3 , ’MT: ' ) 

DO  9  1=1 , NSK 
ISTAT=0 

CALL  MTIN(3, NTAPE, 1000, IA,ISTAT) 

I F ( ISTAT . EQ . 1 )G0  TO  3000 

CONTINUE 

GO  TO  15 

INPUT  FROM  DISC  FILE 

CALL  ASSIGN (3, ’TAN. DAT’) 

DO  11  1=1 , NSK 

READ (3 , END  =  3000, ERR  =  3000  )IA 
CONTINUE 

ASSIGN  OUTPUT  DEVICE 
SET  MTOUT  FOR  D/A 
IF ( MTAPE. LE. 1  )C ALL  ASSIGN (2 ,' MT :' ) 

IF( MTAPE . EQ. 3 )CALL  ASSIGN (2, 'STRESS . DAT  ' ) 

IF ( MTAPE. LE. 1 )CALL  MTSET (2 , MTAPE , 1 , ISTAT) 

BLOCK=0 . 


CALCULATE  PEAK  AMPL  OF  CW  TONES 
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c 

16 

C 

17 
C 

18 
C 


C 


22 


23 

25 

29 


26 

27 

23 


IF (NCW. LE . 0  )G0  TO  16 
AMPCW=SQRT(2./NCW)*10.**(AMPCW/20. ) 

CALCULATE  PEAK  AMPL  OF  FSK  TONES 

IF (NFSK . LE . 0 )GO  TO  17 

AMPFSK=SQRT(2./NFSK)*10.**(AMPFSK/20. ) 

CALCULATE  PEAK  AMPL  OF  PSK  TONES 

IF ( N PSK . LE . 0 )GO  TO  18 

AMPPSK=SQRT(2./NPSK)*10.**(AMPPSK/20. ) 

CALCULATE  PEAK  AMPL  OF  SWEPT  TONE 
AMPSAW=SQRT(2. )» 1 0 .»*( AMPS AW/20 . ) 

CALCULATE  GAUSSIAN  NOISE  SCALE  FACTOR 

Z=10.»AL0G10(SAMP/2. ) 

GAUSS=10.**((GAUSS+Z)/20.  ) 

CALCULATE  NOISE  PULSE  SCALE  FACTOR 
AMPPUL=10.**( (AMPPUL+Z)/20. ) 

CALCULATE  AMPL  SCALE  FACTOR  FOR 
EACH  DELAY  PATH 

CONVERT  DELAY  IN  SEC.  TO  NO.  OF  SAMPLES 

IF (NPATH . EQ. 1 )GO  TO  25 
R  =0  . 

DO  22  1  =  1  ,  NPATH 
R=R+VAR(I ,4)*VAR(I ,4) 

IFCR.LT. 0.001 ) R  =0 . 00 1 
R  =SQRT ( R ) 

DO  23  1=1 , NPATH 

VAR  ( I ,4)=VAR(I ,  4  )  /  R 

CONTINUE 

DO  29  1=1 , NPATH 

VAR(I,  3)  =  VAR(I ,3)*SAMP 

IF  FADING  CHANNEL  READ  IN  PATH  DATA 
PATH  DATA  GENERATED  AT  0.01*SAMP 
STORE  OLD  AND  NEW  VALUES 

IF ( NFADE . EQ . 0 )G0  TO  28 
CALL  ASSIGN ( 4 , ' PATH . DAT  '  ) 

READ (4 , END =3 000, ERR =3000) PATH 
NDY  =0 
DELTA=0 . 

DO  27  J=1 ,2 
DO  26  1=1 , NPATH 
NDY=NDY+1 
ROLDCI )=RNEW( I ) 

GOLD (I  )=GNEW( I  ) 

RNEW(I )= REAL (PATH (NDY) ) 

GNEWCI )=AIMAG( PATH (NDY) ) 

CONTINUE 

CONTINUE 

CONTINUE 


20 


C 

C 


MAIN  LOOP 

CONTINUE 
FRAME=FRAME+1 . 0 

IF  NO  INPUT  , BYPASS  ROUTINES 
WHICH  PROCESS  INPUT  DATA 
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30 

40 


c 


50 


IF(NTAPE.LT.4)G0  TO  40 
DO  30  I=1,NBD 
A ( I ) =0  . 


GO  TO  50 
CONTINUE 


READ  INPUT  DATA 


ISTATiO 
CALL  INPUTX 

IF(ISTAT.EQ. 1)G0  TO  3000 

TRANSMITTER  FILTER 


CALL  TXFILT 


ALL  PASS  FILTER 


CALL  ALPASX 
KSTAT  =  0 


COMBINE  OUTPUTS  FROM  N  PATHS 

CALL  DELAYX 

IF(KSTAT. EQ. 1 )G0  TO  3000 


CALL  FREQX 
CONTINUE 

CALL  INTERF 

CALL  PULSE 

CALL  NOISEX 

CALL  RXFILT 

CALL  AGC 

STORE 

CALL  STOREX 
IF(FRAME.LT. FRAMES  )G0 


FREQ  TRANSLATE  COMPOSITE  SIGNAL 
DOPPLER  FREQ  CHANGES  DUE  TO 
SCENARIO  MADE  IN  SUBROUTINE  DOPVAR 
WHICH  IS  CALLED  FROM  FREQX 


ADD  NARROWBAND  TONAL  INTERFERENCE 
ADD  NOISE  PULSE 
ADD  GAUSSIAN  NOISE 
RECEIVER  FILTER 
RECEIVER  GAIN  CONTROL 
OUTPUT 

TEST  FRAME  COUNT 
TO  20 


LOOP  COMPLETED 

3000  CONTINUE 

WRITE (5 , 300 1 )TOTAL 

3001  FORMAT( IX, 'TOTAL  NO.  OF  SAMPLES:  '  ,F9 . 0 ) 

IF (TOTAL. LE . 1 . )G0  TO  3100 

RMS=SQRT ( RMS/TOTAL ) 

V=RMS*5./2047. 

DB  M  =  1  0  .  *  A  LOG 10(V*V/0.6) 

DB:10.*ALOG 10 (PEAK/RMS) 

WRITE(5, 3002)PEAK,DB ,DBM 

3 002  FORMAT ( IX, ' PEAK: ' ,F9 . 2 , '  PEAK/RMS (DB )=' ,F9 . 2, ’  DBM=’,F8.3) 
3100  CONTINUE 

IF( MTAPE . GT . 1 )G0  TO  3200 

CALL  MTOUT (2 , MTAPE , 1000 , 13 , ISTAT ) 

CALL  MTSET(2, MTAPE, 0, ISTAT) 

CALL  MTOUT (2 , MTAPE ,0> IB , ISTAT) 
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CALL  MTOUT (2 , MTAPE , -1 , IB , ISTAT  ) 

GO  TO  3300 
3200  WRITE (2  )IB 
ENDFILE  2 
3300  CONTINUE 

IF(NTAPE . LE . 1 )CALL  MTIN ( 3 , NTAPE , -1 , IA , ISTAT ) 

IF (NTAPE . EQ . 3 )CALL  CL0SE(3) 

IF ( NF ADE . GT . 0 )C ALL  CLOSE (4  ) 

END 

C  INPUT  ONE  FRAME  OF  ANALOG  SAMPLES 

C 

SUBROUTINE  INPUTX 

COMMON /BLK1 /VARC10, 10) ,A( 100) ,Q(100) 

C0MM0N/BLK8/IA( 1000) , ISTAT , NDX , BLOCK 

READ  IN  NBD  SAMPLES 
CONVERT  TO  FLOATING  POINT 
NORMALIZE  TO  ZERO  DBM( V=0 . 77459667 ) 
E=V*2047/5  =317.12 

NB  ^ =VAR (2  ,  1  ) 

NTAPE  =VAR ( 3 , 2  ) 

DO  10  1  =  1  ,  NBD 
NDX=NDX+ 1 

A( I )=FL0AT(IA(NDX)/16  )  / 3 1 7 .  12 
IFCNDX.LT. 1000  )G0  TO  10 
NDX=Q 
ISTAT=0 

IF ( NTAPE . LE . 1 )C ALL  MTIN (3 , NTAPE , 1000 , IA , ISTAT) 

CHECK  FOR  SUCCESSFUL  READ 
IF ( I STAT . EQ . 0 )GOTO  4 

CHECK  FOR  EOF  ON  TAPE 
IFCISTAT. EQ. 1 )GOTO  30 

CHECK  FOR  END  OF  TAPE 
IFCISTAT. NE.2)G0T0  3 
WRITE (5,22) 

22  FORMATC '  END  OF  TAPE  ' ) 

GO  TO  30 

C  TAPE  READ  ERROR 

3  WRITE (5 , 33 ) ISTAT 

33  FORMATC’  TAPE  READ  ERROR  ISTAT=:’,I5) 

4  IF(NTAPE.EQ.3)HEAD(3,END=30,ERR=30)IA 
BL0CK=BL0CK+1 .0 

C  CHECK  BLOCKER  CODE 

K  =  3 

DO  8  J  =  1 ,  1000 
KB  =  1 A C  J  )  .AND.  15 
IF ( K . EQ . KB )G0  TO  7 
WRITE (5 , 5 ) BLOCK 

5  FORMATC IX, ’3L0CKER  CODE  ERROR  IN  BLOCK’, F8.0) 

GO  TO  10 
7  K=K+4 

IFCK  .GT.  15  )K  =  3 
3  CONTINUE 
10  CONTINUE 
RETURN 
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30  WRITE (5,35)  BLOCK 

35  F0RMAT( IX, ’EOF  OR  ERROR  ON  ANALOG  INPUT  IN  BLOCK  '  ,F8 . 0 ,/ ) 
I STAT= 1 
RETURN 
END 

SUBROUTINE  MTIN (LUN , IDR IVE , ICNT, IARRZ V , ISTAT ) 

MTIN  MAG-TAPE  DRIVER 
D.  TATE  OR  I , INC  JAN  79 

THIS  SUBROUTINE  IS  USED  TO  READ  FROM  MAG  TAPE 
(NON-ANSI).  VARIABLES  ARE  DEFINED  IN  I/O  DRIVERS 
REFERENCE  MANUAL,  SUBROUTINES  ARE  DEFINED  IN 
EXECUTIVE  REFERENCE  MANUAL. 

PARAMETERS: 

LUN  =  LOGICAL  UNIT  NUMBER  (POSITION  IN  LOGICAL 
UNIT  TABLE) 

IDRIVE  =  DRIVE  NUMBER  (0  OR  1) 

ICNT  <  0  =>  REWIND  REQUEST 

>  0  =>  NUMBER  IF  INTEGER  WORDS  IN  RECORD 
IARRAY  =  ARRAY  IN  WHICH  DATA  FROM  TAPE  WILL  BE  STORED 
ISTAT  =  STATUS  RETURNED  TO  CALLING  PROGRAM 
=  0  =>  O.K. 

=  1  =  >  END  OF  FILE  DETECTED 
=  2  =>  END  OF  TAPE  DETECTED 
=  3  =>  DATA  OVERRUN  (BUFFER  SIZE  TOO  SMALL) 

<  0  =>  UNSPECIFIED  QIO$  ERROR 

DIMENSION  I3B ( 2 ) ,IPRL(o) , IARRAY (1  ) 

DATA  ISSUC/ 1/ ,  IOATT/"1  400/ , IEDAA/-8/I0RWD/"2400/ 

DATA  IORLB/"1000/,IEDAO/-1 3/  ,  IEE0F/-10/, IEEOT/-62/ 

C  ASSIGN  LUN  TO  MAG  TAPE  UNIT 

CALL  ASNLUN (LUN , 'MT' , IDRIVE , ISTAT) 

IF(ISTAT.NE. ISSUC)  GOTO  99 
C  CLEAR  PARAMETER  LIST 
DO  10  1=1,6 
10  IPR L ( I  )  =  0 

C  REQUEST  EXCLUSIVE  USE  OF  TAPE  DRIVE 

CALL  QIO(IOATT, LUN, 1 , 0 , ISB , IPRL , ISTAT ) 

C  SEE  IF  QIO  DIRECTIVE  ACCEPTED 

IF(ISTAT.NE. ISSUC)  GOTO  99 

C  QIO  DIRECTIVE  ACCEPTED,  NOW  WAIT  FOR  I/O  TO  FINISH 
CALL  WAITFR ( 1 , ISTAT) 

C  CHECK  WAITFR  DIRECTIVE  ACCEPTANCE 
IF(ISTAT. NE. ISSUC)  GOTO  99 
C  DIRECTIVE  ACCEPTED,  NOW  CHECK  I/O  STATUS 
IF(ISB( 1  )  .EQ. ISSUC)  GOTO  20 
C  IF  DEVICE  ALREADY  ATTACHED,  CONTINUE 
C  SIGN  EXTEND  FOR  ERROR  CODES 
ISTS  =  ISB( 1  )  —  2  5  6 
IF(ISTS. NE. IEDAA)  GOTO  99 
C  CHECK  ICNT  FOR  FUNCTION 
20  IF (ICNT. GE . 0 )  GOTO  30 

C  REWIND  REQUEST 

CALL  QIO (IORWD, LUN ,2,0, ISB, IPRL, ISTAT) 

GOTO  80 

C  NORMAL  READ  RECORD  REQUEST 
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CALL  GETADR(IPRL(1 ),IARRAY(1 )) 

IPRL (2 )=ICNT*2 

CALL  QIO(IORLB,LUN,2,0, ISB.IPRL, ISTAT) 

C  CHECK  DIRECTIVE  ACCEPTANCE  FOR  ALL  FUNCTIONS 
30  IF(ISTAT.NE. ISSUC)  GOTO  99 

CALL  WAITFR (2 , ISTAT  ) 

IFUSTAT.NE.  ISSUC)  GOTO  99 
C  FUNCTION  WAS  ACCEPTTED,  NOW  CHECK  FINAL  STATUS 
ISTAT=256 

IF (I  SBC  1  )  .EQ. ISSUC)  ICNT=ISB(2  )/2 
IFCISBC 1  )  .EQ. ISSUC)  GOTO  99 
C  SIGN  EXTEND  STATUS  FOR  ERROR  CODES 
I3TS=ISB( 1  )  —  2  5  6 
IF ( I  STS . EQ . IEEOF )  ISTAT=257 
IFCISTS.EQ. IEEOT)  ISTAT=258 
IF ( I  STS . EQ . IEDAO)  ISTAT=259 
C  NOT  EOF,  EOT  OR  DAO,  BUT  STILL  AN  ERROR 
IF ( I  ST AT . EQ . 256  )  ISTAT  =  ISB( T ) 

99  ISTAT=ISTAT-256 

RETURN 
END 

C  ALLPASS  FILTER 

C 

SUBROUTINE  ALPASX 

C  ALL  PASS  FILTER 

C  OUTPUT  OF  PHASE  SHIFTER  1  IS  IN  ARRAY  Q 

C  ii  ii  it  ii  2  it  "  "  A 

COMMON /BLK1 /VAR ( 10 , 1 0 ) , A ( 1 00  ) , Q ( 1 00  ) 

COMMON /BLK 9 /SAVIN (3, 2) ,SAV0UT(3 ,2) 

DIMENSION  XI ( 10  1 ,2)  ,X2( 100)  ,HC(3 , 2) 

DATA  HC /-. 60746 U7, .3229095, .8083133, .9418089, .61 321 04, -.0730881/ 
NBD=VAR ( 2 , 1 ) 

DO  15  Iil,  NBD 
15  X2(I)sA(I) 

C  FILTER  NBD  POINTS 
C  LA=1 ,  Q  PHASE 
C  LA  =  2 ,  I  PHASE 

NBD1 =  N  B  C  + 1 
DO  10  LA  =  1  ,2 
DO  12  1=1 , NBD 
12  XI  (1  +  1  ,  1  )  =  X2d  ) 

KIN=2 

DO  30  I  SEC  =  1  ,  3 

KOUT  =KIN 

KIN=KIN+1 

IF (KIN . GT . 2  )KIN  =  1 

XI (1 , KIN )=SAVIN(ISEC,LA  ) 

XI ( 1 , KOUT)=SAVOUT(ISEC,LA  ) 

DO  20  I  =2  ,  NB  D  1 

20  XI (I ,KOUT)=HC( I  SEC, LA)* (XI (I , KIN )+X 1 (1-1 , KOUT) ) -XI (1-1 , KIN  ) 

SAVIN (I  SEC, LA )  =  X1 (NBD 1 , KIN ) 

30  SAV0UT(ISEC,LA)=X1(NBD1,K0UT) 

GO  TO  (40, 50  ), LA 
40  DO  45  I =1 , NBD 
45  3 (I )  =  X 1 (I +1 , KOUT  ) 
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10  CONTINUE 
50  DO  200  1=1 , NBD 
200  A  ( I )=X1 (1+1 , KOUT ) 

3000  RETURN 
END 

ADD  GAUSSIAN  NOISE  TO  ANALOG  SAMPLES 
SUBROUTINE  NOISEX 

COMMON /BLK 1 /VAR( 10, 1 0 ) , A ( 1 00 )  , Q ( 1 00 ) 

COMMON /BLK6 /I  1 , 12 
EQUIVALENCE  ( VAR ( 1 ,  2 ) , GAUSS ) 

NBD  =  VAR ( 2 ,  1  ) 

DO  10  1=1 , NBD 
U  =R AN ( I  1 , 12) 

V  =  R AN ( 1 1  , 12) 

U=GAUSS»SQRT(2 . 0*ABS(AL0G(U ) ) ) 

V  =6 . 23  3 1  35  308  *V 
A(I  )  =  A(I  )  +  U*COS(V ) 

10  CONTINUE 
RETURN 
END 

;  FREQUENCY  TRANSLATION  BY  DOPP 

SUBROUTINE  FREQX 

COMMON/BLK1/VAR(lO,  1 0  )  , A ( 1 00  )  , Q ( 1 00 ) 

COMMON /BLK 5/PATH ( 1024) ,WT, TIMEX 
EQUIVALENCE  ( V A  R ( 1 , 1  )  ,SAMP)  ,(VAR(10,9),DOPP) 

COMPLEX  PATH 
NBD=VAR (2,1 ) 

TW0PI=6. 283135303 
TS=TWOPI /SAMP 
DO  50  1  =  1,  NBD 

:  UPDATE  VARIATIONS  IN  DOPPLER 

CALL  DOPVAR 
WT=WT+DOPP*TS 

IF ( WT . GT . TWOPI )WT=WT-TWGPI 
IF ( WT . LT . -TWOPI ) WT=WT+TWOP I 
50  A  ( I ) = A ( I ) *  COS ( WT )+Q ( I ) *S IN ( WT ) 

RETURN 

END 

C  TAP  DELAY  LINE 

SUBROUTINE  DELAYX 

COMMON /3LX1 /VAR ( 10 ,10),A(100),Q(100) 

COMMON /BLK 3 /DELAY (2 , 200 ) 

C0MM0N/BLK5/PATHC 1024) ,WT, TIMEX 

COMMON /BLK 1 3 /DELTA , ROLD( 10  )  , RNEW( 1 0 ) ,GOLD( 1 0 ) ,GNEW( 1 0 )  ,NDY , KSTAT 
COMPLEX  PATH 
NFADE  =VAR (4 , 1 ) 

IF( NFADE . £Q . 0 ) RETURN 
NBD=VAR (2 , 1 ) 

NPATH  =VAR ( 3 , 1 ) 

:  ADD  DATA  TO  DELAY  LINE,  CLEAR  INPUT 

DO  50  1=1, NBD 
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DELAY ( 1 , I+NBD)=A(I ) 

DELAY(2, I+NBD)=Q ( I ) 

A ( I  )  =  0  . 

50  Q  ( I )  =  0 . 

ADD  SAMPLES  FROM  EACH  PATH 
RELATIVE  AMPL  OF  EACH  PATH  IS  IN  VAR(J,4) 
PATH  GAIN  VALUES  SAMPLED  AT  0.01*SAMP 
RELATIVE  DELAY  OF  EACH  PATH  IS  IN  VAR(J,3) 
DELAY  GIVEN  IN  NO.  OF  SAMPLES 

KSTAT=0 

DO  150  1  =  1  ,  NBD 
DELTA  =DELTA+0 .01 
DO  100  J  =  1  , NPATH 
N =NBD-VAR ( J , 3) 

INTERPOLATE  BETWEEN  PATH  VALUES 
R=ROLD(J )+(RNEW( J  )-ROLD(J ) ) ‘DELTA 
G=GOLD ( J )+ (GNEW( J 1 -GOLD (J )) ‘DELTA 
X=DELA Y ( 1 , N  +  I ) ‘R+DELA Y (2 , N  +  I  )‘G 
A(I)=A(I )+X*VAR(J,4) 

X=DELAY(2 , N+I )*R-DELAY( 1 , N+I)*G 
Q ( I )=Q(I )  +  X  ‘  V  A  R  ( J , 4) 
no  CONTINUE 

IF( DELTA . LT . 1 . )G0  TO  150 

READ  IN  NEW  PATH  VALUES 

DELT  A  =0  . 

DO  110  J  =  1 , NPATH 

NDY  =NDY+1 

ROLD ( J  )  =  RNEW( J  ) 

GOLD ( J  )=GNEW( J  ) 

RNEWCJ  )  =  REAL ( PATH (NDY  )  ) 

GNEW(J )=AIMAG( PATH (NDY)) 

IF ( NDY . LT . 1024)G0  TO  1  10 
NDY  =0 

READ (4 , END  =  3 000 , ERR  =  3000  )PATH 
HO  CONTINUE 
150  CONTINUE 

SHIFT  SIGNAL  IN  DELAY  LINE 

DO  200  1=1, NBC 

DELA Y( 1 , I)=DELAY( 1 , I+N3D) 

200  DELAY(2 , I)=DELAY(2, I+N3D) 

RETURN 

3000  WRITE (5 , 4000  ) 

4000  FORMATOX, 'EOF  ON  PATH  DATA',/) 

KSTATs  1 

RETURN 

END 


UPDATE  VARIATION  IN  DOPPLER 
3  ALGOTHRIMS:  STE P , RAM P , AND  FLY3Y 
SUBROUTINE  DOPVAR 

COMMON /BLK1 /VAR ( 1 0 , 10  )  , A( 1 00 )  ,Q( 1 00 ) 

COMMON /BLK 5 /PATH( 1024 ) ,WT, TIMEX 
COM MON/8 LX  7/FLAG,. KF LAG 

EQUIVALENCE  (VAR (1,1)  ,SAMP)  ,  ( VAR ( 1  , 9 ) , START )  ,  ( VAR (2 , 9 ) ,DOPX) 
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EQUIVALENCE  (VAR (3 , 9) , RATE ) , ( VAR ( 4 , 9 ) .DOPLOW) , (VAR (5 ,9) , DOPHGH ) 
EQUIVALENCE  (VAR (6 , 9 ) , TIMEY ) , ( VAR (7 , 9 ) , XDIST  )  ,  ( VAR ( 3 , 9  )  ,  VEL  ) 
EQUIVALENCE  ( V AR (9 , 9 ) , RF ) , ( V A R ( 1 0 , 9  )  ,  DOPP  ) 

COMPLEX  PATH 
C 

KDOPP  =VAR (5 , 1 ) 

IF(KD0PP.EQ.0)RETURN 
GO  TO  (1 ,2, 3), KDOPP 

C  STEP  CHANGE  IN  DOPPLER  AT  TIME  "START" 

C  TEST  IF  STEP  HAS  OCCURRED 

1  IF (KFLAG. EQ.  1  )  RETURN 
TIMEXsTIMEX+1 .0/SAMP 
IF ( TIMEX. LT. START)  RETURN 
DOPP=DOPP+DOPX 
KFLAGsl 
RETURN 

RAMP  CHANGE  IN  DOPPLER  AT  CONSTANT  RATE 
MAX  LIMIT=  DOPLOW, DOPHGH 

2  DOP=DOP+RATE/SAMP 
IF( DOPP . GE . DOPHGH) RATE =-1 . 0*RATE 
IF (DOPP . LE. DOPLOW) RATE=- 1 . 0*RATE 
RETURN 

HIGH  VELOCITY  FLYBY 
BEGIN  TEST  AT  TIME  START 
TIMEY  IS  TIME  PLANE  IS  FROM  MIN  DISTANCE 
XDIST  IS  MIN  DISTANCE 
VEL  IS  MPH 
RF  IS  MHZ 

3  riMEX=TIMEX+ t .O/SAMP 
IF (TIMEX.LT. START) RETURN 
V= (VEL/3600 . ) **2 
X=XDIST*XDIST 
TsTIMEY*TIMEY 
XI =SQRT(X+V*T) 

DOPPsV»RF*TIMEY/X 1 /O . 186 
TIMEYsTIMEY-1 .O/SAMP 
RETURN 
END 


ADD  MULTITONE  INTERFERENCE 
CW  ,  FSK  ,  4-PHASE  DPSK,  SWEPT  TONE 
MAX  OF  40  TONES  IN  FSK  AND  DPSK  MODES 

SUBROUTINE  INTERF 

COMMON /BLK 1 /VAR( 10 , 10),A(100),Q(100) 

C OMM ON /BLK 4/PHASE (40)  , TAU ( 4  )  , FQS AW 
C0MM0N/BLK6/I  1  , 12 
DIMENSION  GR AY ( 4  ) 

EQUIVALENCE  (VAR (1,1) , SAMP ) , ( VAR (2 , 5 )  , F 1 CW)  , ( V A R ( 3 , 5 ) , FXCW ) 
EQUIVALENCE  (VAR (4 , 5 ) , AMPCW)  , ( V A R (2 , 6 ) , F 1MK ) , ( V AR ( 3 , 6  )  , SH IF T  ) 
EQUIVALENCE  (VAR ( 4 , 6 ) , FXFSK ) , ( V AR (5 , 6 ) , AMPFSK ) ,(VAR(2,7),F1PSK) 
EQUIVALENCE  ( V  A R ( 3 , 7 ) , FXPSK ) , ( V AR ( 4 , 7 ) , AMPPSK )  ,  ( V A R ( 1 , 8 ) , F 1 S AW) 
EQUIVALENCE  (VAR (2 , 8 ) ,F2SAW)  , (VAR (3 , 3 )  ,  RSAW)  , (VAR (4 ,8) , AMPS AW) 
DIMENSION  FFSK ( 40  ) 
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DATA  GRAY/1.0,3.0,-1.0,-3.0/ 

KINTF  =VAR (6  ,  1  ) 
IF(KINTF.EQ.0)RETURN 
PI =3  -  14159 
TWOPI =2 . 0#PI 
NBD=VAR (2  ,  1  ) 

MCW=VAR ( 1 ,5) 

NFSK=VAR ( 1 ,6) 

NPSK=VAR ( 1 ,7) 

TS= 1 ./SAMP 
M=KINTF . AND . 1 
I F  (  M  .  EQ .  0  )GO  TO  2 


DO  10  1=1 , NBD 
FREQ=F 1CW-FXCW 
TAU(1  )  =  T  A  U ( 1  )+TS 
IF(TAU(1 ) . GE . 1 . 0 )T AU( 1  )  =  TAU(1  )-1  . 

DO  5  J  =  1 , NOW 
FREQ=FREQ+FXCW 
W=TWOPI *FREQ 

5  A ( I ) =A ( I )+AMPCW*COS ( W*T  AU ( 1  )) 

10  CONTINUE 
C 

C  BINARY  FSK 

2  M =KINTF . AND . 2 

IF ( M . EQ. 0 )G0  TO  3 

"NFSK"  TONE  PAIRS 
FSK  SHIFT="SHIFT"  HZ 
SPACING  BETWEEN  ADJACENT  PAIRS  IS 
"FXFSK"  HZ,  PEAK  AMPL  OF  EACH  TONE  IS 
" AMPF3K"  ,  RANDOM  KEYING 

FREQ=F  1MK-FXFSK 
DO  60  1  =  1  ,  NFSK 
FREQ=FREQ+FXFSK 
FFSKCI  )  =  FREO 
X  =  RAN ( I  1,12) 

IF(X.GT.0.5)FFSK(I)  =  FFSK(I  )  +  SHIFT 
60  CONTINUE 

DO  20  1  =  1  ,  NBD 
FREQ=F 1MK-FXFSK 
TAU(2 )  =  TAU(2  )+TS 
IF(TAU(2).GE.1.)TAU(2)=TAU(2)-1. 

DO  15  J=1 , NFSK 
W  =TWOP I *FFSK ( J  ) 

15  A ( I ) = A ( I )+AMPFSK*COS (W*TAU(2)) 

20  CONTINUE 


CW  INTERFERENCE 

" NCW"  TONES  SPACED  "FXCW"  HZ  APART 
PEAK  AMPL  OF  EACH  TONE  IS  "AMPCW" 


3  M  =  KINTF . AND.  4 
IF ( M . EQ . 0 )G0  TO 


DPSK ,  4PHASE 

"NPSK"  TONES  SPACED  '’FXPSK"  HZ  APART, 
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PEAK  AMPL  OF  EACH  TONE  IS  "AMPPSK" 
RANDOM  KEYING 

DO  70  1=1 , NPSK 
l  =  1 

X  =  RAN ( I  1 , 12) 

IF  (  X  .  GT  .  0 . 5  )L  =2 
X  =  RAN ( I  1,12) 

IF ( X . GT . 0 . 5 )L  =1  +  1 

PHASE(I)=PHASE(I)+GRAY(L)*PI/4.0 
CALL  LIMIT (PHASE ( I ) ) 

CONTINUE 
DO  30  1=1 , NBD 
FREQ=F 1PSK-FXPSK 
TAU(3)=TAU(3 )+TS 

IF(TAU(3).GE.1.0)TAU(3)=TAU(3)-1  .0 

DO  25  J=1 , NPSK 

freq=freq+fxpsk 

W=TWOPI*FREQ 

A( I  )=A( I  )  +  AMPPSK*COS(W*TAU( 3 )  +  PHASE(J ) ) 

CONTINUE 


CONTINUE 
M=KINTF.AND.8 
IF ( M . EQ . 0  )RETURN 


SWEPT  TONE  FQSAW  BETWEEN  F1SAW  AND  F2SAW 
SWEEP  RATE  IS  RSAW 
AMPL  IS  AMPSAW 


DO  40  1=1, NBD 

FQSAW=FQSAW+R SAW/ SAMP 

IF ( FQSAW .GT . F2S AW) FQSAW =F 1SAW 

TAU(4  )=TAU(4 )+TS 

IF(TAU(4).GE. 1 . )TAU(4)=TAU(U)-1 

W=TWOPI*FQSAW 

A ( I  )  =  A ( I  )+ AMPS AW* COS ( W*T  AU ( 4 ) ) 

RETURN 

END 


ADD  NOISE  PULSE  TO  SIGNAL 

SUBROUTINE  PULSE 

COMMON /BLK 1 /VAR( 10, 10),A(100),Q(100) 

C0MM0N/BLK6/I 1 , 12 
COMMON /BLK 7 /FLAG , KF LAG 

EQUIVALENCE  (VAR ( 1  , 1  )  ,  SAMP ) , ( V A R ( 1  ,  10)  ,PROB) 
EQUIVALENCE  ( V A R ( 2  ,  1 0  )  ,  AMPPUL )  ,  ( V  A R ( 3 , 1 0 ) .PULDUR) 
KPULSE=VAR (7 , 1 ) 

IF ( KPULSE . EQ. 0  5RETURN 


DUMMY  SUBROUTINE 
I  PUT  SIGNAL  IS  3  F.  A  L  !*.’  A-RAY  A 

j  t-r*  pi  i  -y  •»  **  »*  M  ”  M 

V-’PP!JL=  RMS  AMPLITUDE  OF  ‘'OISE  PULSE 
PULDUR  =  DURATION  OF  NOISE  PULSE 


o  o  o  o  o  o  o  o  o  o  o 


PROP  =  OROBAB ILITY  OF  OCCURRENCE  OF  NOISE  PULSE 
FLAG  =  CONTROL  VARIABLE  FOR  NOISE  PULSE 
SAMP:  3 AM  PL  I NO  RATE 

NB  D  =  NO.  OF  SAMPLES  PROCESSED  PER  BLOCK 
II  ,12=  RANDOM  NUMBER  GENERATOR  VARIABLES 


RETURN 

END 


FILTER  S IGN AL+NOISE 

SUBROUTINE  RXFILT 

COMMON /BLK 1 /VAR(10 , 10 ) , A(  1  00 ) ,Q ( 1  00 ) 

COMMON /BLK10/RXFIL(200) 

K=VAR (9 , 1 ) 

IF (K . LE . 0 ) RETURN 

IF  K=1,  BASIC  FILTER  IS  .5Z0-.5Z2 
IF  K=2,  ITERATE  FILTER  1  TWO  TIMES 
IF  K= 3,  FILTER  IS  . 6Z0- . 5Z2- . 1Z 4 
IF  K=4, ITERATE  FILTER  3  TWO  TIMES 
IF  K=5 ,  ITERATE  FILTER  3  FOUR  TIMES 

NBD=VAR (2 , 1 ) 

IF(K . GT . 2 )G0  TO  9 

BASIC  FILTER  IS  .5Z0-.5Z 2 
ITERATE  K  TIMES 

DO  8  I  =  1  ,  K 
DO  2  J  =  1 , NBD 
2  R XFIL ( J  +2  )  =A ( J  ) 

DO  4  J=1 , NBD 

4  A( J )= . 5*RXFIL(J+2 . 5*RXFIL(J  ) 

DO  5  J  =  1  ,2 

5  RXFIL(J ) = RXF IL ( J +NBD ) 

8  CONTINUE 
RETURN 

9  K  =K -2 
IF(K  . GT . 2  )K  =4 

BASIC  FILTER  IS  . 6Z 0- . 5Z 2- .  1 Z 4 
ITERATE  FILTER  K  TIMES 

DO  100  1=1 ,K 
DO  10  Jsl  ,  NBD 
10  RXFILCJ+4 )=A( J ) 

DO  20  J  =  1  , NBD 

20  A(J)=0.6*RXFIL(J+4)-0.B*RXFIL(J+2)-0.1*RXFIL(J) 

DO  30  J  =  1  ,  4 

30  RXFIL(J )=RXFIL(J+NBD) 

100  CONTINUE 
RETURN 
END 

C  TX  FILTER 

SUBROUTINE  TXFILT 

COMMON /BLK 1 /VAR ( 1 0 , 10 )  ,  A ( 1 00  )  , Q ( 1 00  ) 

COMMON /BLK2/TXFIK  200) 
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K=VAR  (8 ,  1  ) 

IF(K . EQ. 0)RETURN 

C  BASIC  FILTER  IS  .6Z0-.5Z2-.1Z4 

C  ITERATE  K  TIMES 

NBD=VAR (2 , 1 ) 

DO  ICO  1  =  1  , K 
DO  10  J=1 , NBD 
10  TXFIL ( J+4 )=A ( J ) 

DO  20  J=1 , NBD 

20  A(J)=0.6*TXFIL(J+4)-0.5*TXFIL(J+2)-0. 1*TXFIL(J  ) 

DO  30  J=1 , 4 

30  TXFIL(J )=TXFIL(J+NBD) 

100  CONTINUE 
RETURN 
END 

SUBROUTINE  MTOUT (LUN , IDR IVE , IC NT , IARRAY , ISTAT ) 

MTOUT  MAG-TAPE  DRIVER 
D.  TATE  OR  I , INC  JAN  79 

THIS  SUBROUTINE  IS  USED  TO  WRITE  TO  MAG  TAPE 
(NON-ANSI).  VARIABLES  ARE  DEFINED  IN  I/O  DRIVERS 
REFERENCE  MANUAL,  SUBROUTINES  ARE  DEFINED  IN 
EXECUTIVE  REFERENCE  MANUAL. 

PARAMETERS: 

LUN  =  LOGICAL  UNIT  NUMBER  (POSITION  IN  LOGICAL 
UNIT  TABLE) 

IDRIVE  =  DRIVE  NUMBER  (0  OR  1) 

IC  NT  <  0  =>  REWIND  REQUEST 

=  0  =  >  WRITE  EOF  REQUEST 
>  0  =>  NUMBER  IF  INTEGER  WORDS  IN  RECORD 
IARRAY  =  ARRAY  TO  TRANSFER  TO  TAPE 
ISTAT  =  STATUS  RETURNED  TO  CALLING  PROGRAM 
=  0  =  >  O.K. 

=  1  =  >  WRITE  PROTECTED  (NO  RING) 

=  2  =>  END  OF  TAPE  DETECTED 
<  0  =>  UNSPECIFIED  QIO$  ERROR 


DIMENSION  IS8(2 ) ,IPRL(6 )  ,IARRAY( 1  ) 

DATA  ISSUC/ 1 /, IOATT/" 1400/, IEDAA/-8/I0RWD/”2400/ 
DATA  IOEOF/"3DOO/ , IOWLB/"400/ , IEWLK/- 1 2/ , IEEOT/-62/ 
C  ASSIGN  LUN  TO  MAG  TAPE  UNIT 

CALL  ASNLUN (LUN , ' MT ' , I  DR  I VE , ISTAT) 

IFdSTAT.NE.  ISSUC)  GOTO  99 
C  CLEAR  PARAMETER  LIST 
DO  10  1=1,6 
10  IPRL(I )=0 

C  REQUEST  EXCLUSIVE  USE  OF  TAPE  DRIVE 

CALL  QIO( IOATT, LUN, 1 , 0 , ISB , IPR L , ISTAT ) 

C  SEE  IF  QIO  DIRECTIVE  ACCEPTED 

IFdSTAT.NE.  ISSUC)  GOTO  99 

C  QIO  DIRECTIVE  ACCEPTED,  NOW  WAIT  FOR  I/O  TO  FINISH 
CALL  WAITFR ( 1 .ISTAT) 

C  CHECK  WAITFR  DIRECTIVE  ACCEPTANCE 
IFdSTAT.NE.  ISSUC)  GOTO  99 
C  DIRECTIVE  ACCEPTED,  NOW  CHECK  I/O  STATUS 
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T  ST-T<!B  /  t  ^ 

IF( ISB( 1  )  .EQ. ISSUC)  GOTO  20 
IF  DEVICE  ALREADY  ATTACHED,  CONTINUE 
SIGN  EXTEND  FOR  ERROR  CODES 
ISTS=ISTAT-256 
IF(ISTS.NE. IE DA A)  GOTO  99 
C  CHECK  ICNT  FOR  FUNCTION 
20  IF(ICNT.GE.O)  GOTO  30 

C  REWIND  REQUEST 

CALL  QIO (IORWD, LUN , 2, 0 , ISB , IPRL , ISTAT ) 
GOTO  80 

30  IF ( ICNT . GT  .  0  )  GOTO  4C 

C  WRITE  EOF  REQUEST 

CALL  QIO (IOEOF, LUN, 2,0, ISB, IPRL, ISTAT) 
GOTO  80 

C  NORMAL  WRITE  RECORD  REQUEST 
40  CALL  GETADR ( I PRL ( 1 ) , I ARRA Y ( 1 ) ) 

IPRL(2 )=ICNT*2 

CALL  QIO ( IOWLB , LUN, 2,0, ISB, IPRL, ISTAT) 

C  CHECK  DIRECTIVE  ACCEPTANCE  FOR  ALL  FUNCTIONS 
80  IFdSTAT.NE.  ISSUC)  GOTO  99 

C A«L  WAITFRC2, ISTAT) 

IF(ISTAT.NE. ISSUC)  GOTO  99 
C  FUNCTION  WAS  ACCEPTTED ,  NOW  CHECK  FINAL  STATUS 
ISTAT=0 

IFdSB(  1  ).EQ. ISSUC)  GOTO  999 
C  SIGN  EXTEND  STATUS  FOR  ERROR  CODES 
ISTS=ISB( 1  )-256 
IFdSTS.  EQ.IEWLK)  ISTAT  =  257 
IF (I STS . EQ . IEEOT )  IoTAT=258 
C  NOT  LOCKED  OR  EOT,  BUT  STILL  AN  ERROR 
IF (ISTAT . EQ . 0  )  IiTAT=ISB(1) 

C  ERROR  RETURN 

99  ISTAT=ISTAT-256 

C  NORMAL  RETURN 

999  RETURN 

END 

C  LIMIT  THE  ANGLE  P  TO  +-PI 

SUBROUTINE  LIMIT(P) 

IF ( P . LT . 0 . 0 )G0  TO  100 
C  P  IS  POSITIVE,  LIMIT  TO  +-TW0PI 
DO  10  1=1 ,500 

IF ( P . LT . 6 . 2831 853072 )G0  TO  20 
10  P=P-6. 2831 853072 

C  LIMIT  P  +-PI 

20  IF  (  P . GT . 3 •  141592  6536)P  =  P-6.28  31  853072 
GO  TO  300 

C  P  IS  NEGATIVE,  LIMIT  TO  -TWOPI 
100  DO  150  1=1 , 500 

IF ( P . GT . -6 . 28 3 1  853072  )G0  TO  160 
150  P=P+6 .2831 853072 

C  LIMIT  TO  +-PI 

160  IF ( P . LT . -3 . 14 15926536 )P=P+6. 2831853072 
300  RETURN 
END 
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SUBRCUTINE  AGC 

COMMON  /PLK  1  /VAR ( 1 0 ,  1  0  )  ,  A  (  1  00  )  ,C  (  P3  ) 

OCVMCM/BLK 1 1 /AGCF  TL (2 ,  100) 

EQUIVALENCE  ( VAR  (2  ,  ?  ) ,  GA  IM  >  ,  ( VA  R  ( 8  ,  2 )  ,  DG  C  A  Y )  ,  ( V  A  E  (7 , 2  )  ,  AGCVI* 


KAGCsVAR (10, 1  ) 

IF ( KAGO  .  EQ. 0  )  RETURN 
Mir  =  VAR(2,  1  ) 

N  A  OC  =V  A  R  ( 5 , 2  ) 


C  DUMMY  SUBROUTINE 

C  GA IM  =  CURRENT  GAIN  VALUE 

C  MAG Os  ATTACK  CONSTANT 

C  DEC  AY=  TOC  AY  CONSTANT 

0  AGC MI. Ms  CONTROLS  max  AMPLIFICATION 


RETURN 
END 

SUBROUTINE  MTSET (LUN , IDRIVE , ISET, ISTAT) 

MTSET  MAG  TAPE  DRIVER  (SET  CHARACTERISTICS) 

D.  TATE  OR  I, INC  AUG  80 

THIS  SUBROUTINE  SETS  THE  TAPE  CHARACTERISTICS.  IT 
DOES  NOT  ATTACH  THE  UNIT.  VARIABLES  ARE  DEFINED  IN 
THE  I/O  DRIVERS  REFERENCE  MANUAL,  SUBROUTINES  ARE  DEFINED 
IN  THE  EXECUTIVE  REFERENCE  MANUAL. 

PARAMETERS: 

LUN  -  LOGICAL  UNIT  NUMBER 
IDRIVE  -  MAG  TAPE  DRIVE  NUMBER  (0  OR  1) 

ISET  -Os  SET  DRIVE  FOR  NORMAL  REWRITE  CAPABILITY 
1  s  SET  DRIVE  FOR  NO  REWRITF  CAPABILITY 
NOTE:  ISET  SHOULD  3E  SET  TO  '  THE  TAPE  IS  TO  BE 
PLAYED  IN  THE  ANALOG  MODE.  THE  ANALOG 
CONTROLLER  CAN  NOT  TOLERATE  EXTENDED  INTER¬ 
RECORD  GAPS.  MTOUT  WRITE  ERRORS  WILL  RETURN 
WITH  ISTAT  =  -56  (IE.BBE  -  BAD  BLOCK) 

ISTAT  -  ERROR  STATUS  RETURNED  TO  CALLING  PROGRAM 
(0  =  NO  ERROR) 

DIMENSION  I  PR L ( 6  )  ,  ISB(2  ) 

DATA  IPRL/6*0/, IS8/2  *0  / 

DATA  IOSTC/"2500/ , ISSUC/1  / 

C 

C  ASSIGN  LUN  TO  MAG  TAPE  UNIT 

CALL  ASNLUN (LUN , ’MT’ .IDRIVE, ISTAT) 

IFdSTAT.  NE.  ISSUC)  GOTO  99 
C  DETERMINE  REQUESTED  CHARACTERISTIC 
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IF ( I  SET . NE . 0 )  IPRL ( 1 )=128 
C  SET  TAPE  CHARACTERISTICS 

CALL  QIO(IOSTC,LUN, 1 , 0, ISB , IPRL, ISTAT) 

C  SEE  IF  QIO  DIRECTIVE  ACCEPTED 

IF(ISTAT.NE. ISSUC)  GOTO  99 

C  QIO  DIRECTIVE  ACCEPTED,  NOW  WAIT  FOR  I/O  TO  FINISH 
CALL  WAITFR ( 1 , ISTAT) 

C  CHECK  WAITFR  DIRECTIVE  ACCEPTANCE 
IFCISTAT.NE. ISSUC)  GOTO  99 
C  DIRECTIVE  ACCEPTED,  CHECK  I/O  STATUS 
IFCISBO  ).NE.  ISSUC)  GOTO  99 
ISTAT=0 

99  RETURN 

END 

C 

SUBROUTINE  STOREX 

C  STORE  DATA  ON  MT:  OR  DISK 

COMMON/ 8 LK 1 /VAR ( 1 0 , 10 ) , A( 1 00 ) , Q ( 1 00 ) 

COMMON/BLK12/IB( 1000) , JSTAT , NDZ , KCODE , RMS , PEAK , TOTAL 
NBD=VAR (2 , 1  ) 

MTAPE  =V AR ( 4 , 2  ) 

DO  100  1=1 , NBD 

REMOVE  NORMALIZATION 
ACCUMULATE  RMS,  PEAK  VALUE,  TOTAL  NO. 
SAMPLES 

A(I  )  =  A ( I )  *3  1  7  .  12 
RMS=RMS+A( I )*A( I ) 

IF(A(I).GT.P£AK)PEAK=A(I) 

TOTALsTOTAL+1 . 

NDZ  =NDZ  +  1 

C  CONVERT  TO  INTEGER,  ADD  BLOCKER  CODE 

I5(NDZ)=IFIX(A(I))/16+KBC0DE 
KBC0DE=KBC0DE+4 
IF ( KBCODE . GT . 15)KBCODE=3 
IF( NDZ . LT . 1 000 )G0  TO  100 
NDZ  =0 

IF ( MTAPE . GT . 1 )WRITE(2)IB 

IF ( MTAPE. LE. 1 )C ALL  MTOUT ( 2 , MTAPE , 1000 , IB , JSTAT) 

C  CLEAR  INTEGER  ARRAY 

DO  50  Jsl  ,  1000 
50  IB ( J  )  =  0 
100  CONTINUE 
RETURN 
END 
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APPENDIX  B 


PROGRAM  VARGEN 

PROGRAM  TO  GENERATE  TABLE  OF  VARIABLES  FOR  STRESS  PROGRAM 
ALL  ENTRIES  ARE  FROM  KEYBOARD  IN  F9-3  FORMAT  IN  RESPONSE 
TO  STATEMENTS 

ALL  DATA  IS  STORED  IN  IN  ARRAY  VAR ( 1 0 , 10) 

DIMENSION  VAR(10, 10) 

WRITE (5 , 1 ) 

1  FORMATC IX, 'ENTER  SAMP,  BAUD,  PATHS , FADING ',/ ) 

DO  100  1=1,4 
100  READ (5,2  )VAR (1,1) 

2  FORM  AT (F9. 3) 

WRITE ( 5,3) 

3  FORMATC IX, 'ENTER  TYPES-  DOP,  INTERF,  PULSE,  FILTERS (TX , RX  )',/ ) 

DO  200  1=5,9 

200  READ ( 5 i 2 )VAR (I , 1 ) 

WRITE (5,20) 

20  FORMATC IX, 'ENTER  AGO  TYPE,  GAUSSIAN  NOISE  LEVEL',/) 

READC5 , 2 )VAR ( 10 , 1 ) 

READ (5,2 )VAR ( 1 ,2) 

WRITEC5 ,66  ) 

66  F0RMATC1X, 'ENTER  AGO  GAIN,  AGO  ATTACK,  AGC  DECAY,  AGCMAX  '  , / ) 
READ(5,2)VAR(2,2) 

READC5 , 2 )VAR (5,2) 

READC5 , 2)VAR (6 ,2) 

READ(5 , 2 )VAR (7,2) 

WRITE (5 , 67  ) 

67  F0RMATC1X, 'ENTER  INPUT,  OUTPUT  DEVICE',/) 

READ (5 , 2 ) VAR ( 3 , 2 ) 

READ ( 5 , 2  )VAR ( 4 , 2 ) 

N  =VAR ( 3 , 1 ) 

WRITE (5 , 4  )N 

4  FORMATC IX, 'ENTER  DELAYTIMES  FOR  ',12,'  PATHS',/) 

DO  300  1=1 , N 

300  READC5 , 2 )VAR (1,3) 

WRITE (5 , 6 )N 

6  F0RMATC1X, 'ENTER  DELAY  AMP  FOR  ’,12,'  PATHS',/) 

DO  500  1  =  1  ,  N 

500  READ (5 , 2 )VAR (I , 4 ) 

K=VAR (6 , 1 ) 

M=K . AND .  1 

IF ( M . EQ . 0 )G0  TO  22 
WRITE (5,7) 

7  F0RMATC1X, 'ENTER  CW  NTONE , LOWEST ,  SPACING,  AMPL ’  ,/ ) 

DO  600  1=1 ,4 

600  READ ( 5 , 2 )VAR (I , 5 ) 

22  M=K .AND .  2 

IF ( M . EQ . 0 )G0  TO  24 
WRITE (5,8) 

8  F0RMATC1X, 'ENTER  FSK  NTONE , F 1MARK , SHIFT , SPACING , AMPL ’,/ ) 


B-l 


DO  1600  1=1,5 
1600  READ ( 5 , 2 ) VAR ( I , 6 ) 

24  M=K .AND . 4 

IF ( M . EQ . 0 )G0  TO  26 
WRITE ( 5 , 10) 

10  F0RMAT(1X, 'ENTER  DPSK  NTONE , F 1 , SPACING , AMPL  ',/ ) 

DO  700  1=1,4 

700  READ(5,2)VAR(I ,7) 

26  M=K .AND . 8 

IF ( M . EQ . 0 )G0  TO  30 
WRITE ( 5,27) 

27  FORMATdX,  'ENTER  SWEPT  TONE-  LOW  ,  HIGH  ,  RATE  ,  AM  PL  '  ,  /  ) 

DO  23  1=1,4 

28  READ(5,2)VAR(I  ,8) 

30  iC=VAR(5,1) 

IF ( K . EQ . 0 )G0  TO  40 

WRITE (5 , 1 1 ) 

11  FORMATdX, 'ENTER  START  TIME  FOR  DOPPLER  CHANGE',/) 
READ(5,2)VAR(1  ,9) 

IFOC.GT.  1  )G0  TO  35 
WRITE (5 , 12) 

12  FORMATdX, ’ENTER  DOPPLER  STEP  SIZE',/) 

READC5 , 2 )VAR (2 , 9 ) 

GO  TO  40 

35  IF ( K  .  GT . 2  )G0  TO  38 
WRITE (5 , 36 ) 

36  FORMATdX, 'ENTER  DOPP  RAMP  RATE',/) 

READC5 , 2 )VAR (3 , 9 ) 

WRITE ( 5,13) 

13  FORMATdX, 'ENTER  DOPP  LOW  AND  HIGH  LIMIT’,/) 

READ(5 , 2 )VAR (4 , 9 ) 

READC5 , 2 )VAR (5 , 9 ) 

GO  TO  40 
38  CONTINUE 

WRITE (5 , 14) 

14  FORMATdX, 'ENTER  FLYBY  START  TIME  AND  MIN  DISTANCE',/) 
READ(5,2)VAR(5,9) 

READ(5,2)VAR(7,9) 

WRITE (5 , 15 ) 

15  FORMATdX,  'ENTER  FLYBY  VELOCITY  AND  RADIO  FREQ’,/) 
READ(5 , 2 ) V  A  R  C  3 , 9 ) 

READ(5,2)VAR(9,9) 

40  WRITE (5 ,41) 

41  FORMATdX,  'ENTER  FREQ  TRANSLATION  ERROR',/) 
READ(5,2)VAR(10,9) 

IF( VAR (7 , 1  )  .EQ.O. )G0  TO  50 
WRITE (5 , 16 ) 

16  FORMATdX, 'ENTER  PULSE  PROB  ,  AMPL,  DURATION ’,/ ) 

DO  17  1=1,3 

17  READ(5,2)VAR(I, 10) 

50  CONTINUE 

CALL  ASSIGN(2, 'STRVAR.DAT') 

WRITE (2 )VAR 
ENDFILE  2 

CALL  CLOSE (2 ) 

END 
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APPENDIX  C 
PROGRAM  EDVAR 

EDIT  TABLE  OF  VARIABLES  FOR  STRESS  PROGRAM 

USER  SPECIFIES  LOCATION  TO  BE  CHANGED  AND  NEW  VALUE 
LOCATION  IDENTIFIED  BY  FORMAT  12, IX, 12 
TERMINATE  BY  ENTERING  ZERO 
NEW  VALUE  ENTERED  BY  FORMAT  F9.3 
PROGRAM  IDENTIFIES  OLD  VALUE 

DATA  STORED  IN  DISK  FILE  STRVAR.DAT 

DIMENSION  VAR(10, 10) 

CALL  ASSIGN(2, 'STRVAR.DAT') 

READ (2 ,END=3000, ERR =3000) VAR 
CALL  C LOSE ( 2 ) 

10  WRITE (5,1) 

1  FORMATOX,  'ENTER  LOCATION',/) 

READ ( 5 , 2 ) I , J 

2  F0RMAT(I2, IX,  12) 

IF ( I . EQ . 0 )G0  TO  30 
WRITE (5 , 3  ) V AR ( I , J  ) 

3  FORMAT ( 1 X , F9 . 3  ) 

WRITE (5,4) 

4  FORMAT( IX, 'ENTER  CHANGE  ’  ,/) 

READ(5 , 5  )VAR (I , J  ) 

5  FORMAT(F9.3) 

GO  TO  10 

30  CALL  ASSIGN(2, 'STRVAR.DAT' ) 

WRITE(2)VAR 
ENDFILE  2 
CALL  CLOSE ( 2  ) 

3000  CONTINUE 
END 


C-l 


ooooooociooooooooooooo 


APPENDIX  D 


PROGRAM  PATHM 

GENERATE  A  FILE  OF  COMPLEX  VALUES  EQUAL  TO  THE  SUM 
OF  NTONE  SINUSOIDS  OF  DIFFERENT  FREQUENCY,  REPRESENTING  THE 
COMPLEX  GAIN  OVER  ONE  PATH 

MAXIMUM  OF  10  PATHS,  MAXIMUM  OF  30  TONES/PATH 
DO  THIS  FOR  M  PATHS,  INTERLEAVE  RESULTS 
NTONE  FREQUENCIES  FOR  EACH  PATH  ARE  RANDOMLY  SELECTED 
WITHIN  NTONE  INCREMENTS  OVER  TOTAL  SPREAD 
STORE  IN  BLOCKS  OF  1024  COMPLEX  VALUES 

INPUTS  AT  RUN  TIME  ARE: 

MODEM  SAMPLING  RATE 
NO.  OF  PATHS 

DOPPLER  SPREAD  AND  MEAN  FREQUENCY  OFFSET  FOR  EACH  PATH 
NO.  OF  TONES  USED  TO  GENERATE  GAIN  VALUES 
NO.  OF  SECONDS  OF  DATA 

OUTPUT  STORED  IN  DISK  FILE  PATH . DAT 

DIMENSION  WT(10, 30) ,F( 10,30) 

DIMENSION  SPR(IO) , DOP ( 1 0  ) 

COMPLEX  W( 1024) 

WRITE ( 5  ,  1  ) 

1  FORMATC IX, 'ENTER  MODEM  SAMPLING  FREQ’,/) 

READ(5 , 2 )SAMP 

C  TAP  GAIN  VALUES  GENERATED  AT  0.01*SAMP 

SAMP=SAMP/1 00 . 

2  F  ORM AT ( F  9 • 2 ) 

WRITE (5,3) 

3  FORMATOX, ’ENTER  NO.  OF  TONES/PATH',/) 

READ ( 5 , 4  INTONE 

4  F0RMATCI5 ) 

WRITE (5 , 5  ) 

5  FORMATC IX, 'ENTER  NO.  OF  PATHS,  AND  SPREAD  AND  DOP  FOR  EACH’,/) 
READ(5 , 4  )NPATH 

DO  10  1=1 , NPATH 
WRITE (5,6) 

6  FORMATOX,  'ENTER  SPREAD',/) 

READ ( 5 , 7  )SPR ( I ) 

WRITE (5,8) 

8  FORMATOX, 'ENTER  DOPPLER’,/) 

READ ( 5 , 7  )DOP ( I ) 

7  FORMAT(F9.6) 

10  CONTINUE 

WRITE (5,20) 

20  FORMATC IX, 'ENTER  NO . OF  SECONDS  OF  DATA  TO  GENERATE’,/) 

READC5 , 2)SEC 

CONVERT  SECOND  TO  BLOCKS  OF  1024 
NBLK=1+SEC*SAMP*NPATH/1024. 

11=0 


30 


32 

35 

40 

4  7 

45 

50 


100 


400 


1000 


12=0 

TWOPI *6. 28318 
DO  30  1  =  1  ,  100 
X  =  R  A  N  ( I  1 , 12) 

GENERATE  NTONE  RANDOM  FREQUENCIES  FOR  EACH  PATH 

DO  40  1=1 , NPATH 
FX=2.*SPR(I)/NT0NE 
FY=-SPR (I ) 

DO  35  J=1 , NTONE 
X  =  RAN (I  1,12) 

F(I,J)=FY+X*FX+DOPII) 

X  =  R  A  N ( 1 1 ,12) 

WT(I , J)=X*TWOPI 
WRITE (5 ,32)I,J,F(I,J) 

FORMAT(1X,2I5,2X,F8.4) 

F ( I  ,  J)  =  F(I , J)*TWOPI/SAMP 

FY=F Y+F X 

CONTINUE 

DO  50  1=1 , NPATH 

DO  45  J=1 , NTONE 

IFIABSIFU,  J))  .LT.TWOPDGO  TO  45 

IF(F(I,J).GE.TWOPI)F(I,J)=F(I,J )-TWOPI 

IF(F(I,J).LE.  -TWOPI  )F(I,J)  =  F(I,J  )+TWOPI 

GO  TO  41 

CONTINUE 

CONTINUE 

CALL  ASSIGN ( 3 , ' PATH . DAT ' ) 

NDX=0 

JBLK=0 

SQNT=SQRT( FLOAT (NTONE)) 

START  OF  MAIN  LOOP 
TO  GENERATE  GAIN  VALUES, 

FOR  EACH  PATH  GENERATE  AN  I  AND  Q  VALUE 

CONTINUE 

DO  1000  1  =  1  , NPATH 
A  =0  . 

Q=0 . 

DO  400  Jsl  ,  NTONE 

WT  ( I  ,  J  )  ='WT  ( I  ,  J  )  +F  ( I  ,  J  ) 

IF( WT( I , J ) .GE . TWOPI ) WT(I , J )=WT(I , J ) -TWOPI 
IF ( WT ( I , J ) . LE . -TWOPI )WT(I,J)=WT(I,J )+TWOPI 
A  =A  +COS ( WT ( I , J)) 

Q=Q+3IN(WT(I , J ) ) 

NORMALIZE 

A=A/SQNT 

Q=Q/SQNT 

STORE 

NDX=NDX+1 

W(NDX)=CMPLX(A,Q) 

IF ( NDX . LT . 1024 )G0  TO  1000 
N  D  X  =  0 
WRITE ( 3  )W 

INCREMENT  3 LOCK  COUNT 

JBLK= JBLK+ 1 
CONTINUE 


TEST  BLOCK  COUNT 

IF( JBLK.LT.NBLK)GO  TO  100 
ENDFILE  3 
CALL  CLOSE (3) 

END 


D-2 


